Universal Four-Boson States in Ultracold Molecular Gases: Resonant Effects in 
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We study the manifestations of universal four-body physics in ultracold dimer-dimer collisions. We 
show that resonant features associated with three-body Efimov physics and dimer-dimer scattering 
lengths are universally related. The emergence of universal four-boson states allows for the tunability 
of the dimer-dimer interaction, thus enabling the future study of ultracold molecular gases with both 
attractive and repulsive interactions. Moreover, our study of the interconversion between dimers 
and Efimov trimers shows that B2 + B2 ^ B3 + B rearrangement reactions can provide an efficient 
trimer formation mechanism. Our analysis of the temperature dependence of this reaction provides 
an interpretation of the available experimental data and sheds light on the possible experimental 
realization of rearrangement processes in ultracold gases. 

PACS numbers: 31.15.xj, 21.45.-v, 34.50.-s, 34.50.Cx, 67.85.-d 
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In both nuclear and atomic physics, the simplest ex- 
amples of quantum halo states [J Q are weakly bound 
dimers, with large radii extending well into classically for- 
bidden regions. A remarkable consequence of their large 
size is that quantum-halo dimers obey universal scaling 
laws, i.e., many of their properties are independent of the 
details of their short-range interaction, typically char- 
acterized by a short length scale rg. The notion of a 
halo state extends in a nontrivial way to more complex 
few-body systems, of which three-body Efimov states Q 
are the most prominent example. Recently, the notion 
of universality has been extended to four-boson systems 
0, S, 0| ^"^d a new universal picture has emerged. For 
each Efimov trimer precisely two four-boson states have 
recently been shown to exist, whose energies are univer- 
sally related to the Efimov trimer energy 0, @1- Thus, 
the four-body system inherits many of the characteris- 
tics of three-body Efimov physics, as the geometric scal- 
ing of energies and length scales More recently, a 
first experimental evidence of such universal four-bod. 
states has been found in an ultracold gas of Cs atoms 
through the observation of resonant losses due to four- 
body recombination satisfying the universal predictions 
of Ref. 0. These findings, along with the rapid exper- 
imental advances in controlling few-body correlations in 
ultracold quantum gases [1, 0, [l^, [lH , magnify the im- 
portance of four-body scattering processes. They offer 
a path for the observation of universal physics and can 
potentially enrich the range of experimentally accessible 
phenomena. 

In this Letter, we explore dimer-dimer scattering pro- 
cesses and analyze the consequences of the recently- 
discovered universal properties of four-boson systems 
0, @] ■ Our interest concentrates on such processes near a 
Feshbach resonance Q where, through application of an 
external magnetic field, the s-wave two-body scattering 
length a can be tuned from —00 to +00. Much of the re- 



cent progress in understanding dimer-dimer correlations 
in ultracold physics has been devoted to dimers formed 
from fcrmionic atoms |12| . For bosons, the relatively 



simple fermionic few-body physics is replaced by a much 
more complex structure involving multiple few-body halo 
states that can strongly affect dimer-dimer correlations 
and the collisional behavior of entire gas. Here we advo- 
cate that ultracold molecular gases are perhaps the best 
candidates for exploring universal four-body physics at 
its full complexity. In fact, the search has already begun. 
In Ref. [lO|, Ferlaino et al. formed an ultracold sam- 
ple of Cs2 dimers and found by varying a that the loss 
coefficient exhibits a pronounced minimum which allows 
for longer lifetimes. The temperature dependence of the 
loss rate was also found to display an intriguing behavior, 
deviating from the expected Wigncr threshold law. 

In the present exploration of four-boson universal 
physics we find that, in contrast to the fcrmionic case 
where dimers interact repulsively [l^ , bosonic dimers dis- 
play resonance effects caused by four-body states 0, Q . 
Near such resonances, the dimer-dimer scattering length 
add can vary from — c» to -l-oo, thereby resulting effec- 
tively in either attractive or repulsive interactions. Our 
studies of dimer-trimer conversion, similar in spirit with 
the one in Ref. , further indicate that a molecular gas 
can be used to form one of Efimov trimers. We find that 
when the atom-trimer collision threshold becomes nearly 
degenerate with the dimer-dimer threshold, most of the 
molecular losses are due to i?2 + -B2 ^ B3 + B rear- 
rangement scattering. Because of the vanishingly small 
kinetic energy released in such process, it allows trimers 
to remain trapped. Here we show that the temperature 
dependence of this rearrangement reaction rate can par- 
tially explain some recent experimental observations by 
Ref. which indicate the existence of a trimer state 
just above the dimer-dimer collision threshold. 

We study dimer-dimer processes using the adiabatic 
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hypersphcrical representation for the four-body problem 
jl4| , which offers a simple and conceptually clear descrip- 
tion of the bound and scattering properties. Here, the 
Schrodinger equation reduces to a simple system of ordi- 
nary differential equations given by, 



E 



(1) 



where the hyperradius R describes the overall size of the 
system, and v va a, collective index that represents all 
quantum numbers necessary to label each channel. In 
Eq. Il]), /i is the four-body reduced mass. E is the total 
energy, and F^ is the hyperradial wave function. W^y' (R) 
are the nonadiabatic couplings that drive inelastic tran- 
sitions between channels and W,y(R) arc effective poten- 
tials that support bound and quasi-bound states dictat- 
ing many of the scattering properties of the system. 
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FIG. 1: (color online), (a) Spectrum of the four-boson system 
from our numerical calculations and (b) a schematic rep- 
resentation of the region important for dimer-dimer collisions 
(see text). The four-boson states are represented by black 
solid lines, while the collision thresholds are represented by 
solid red lines for the atom-atom-dimer and dimer-dimer col- 
lisions, and by green dashed lines for atom-trimer collisions. 

Figure [T^a) shows the spectrum for the two-, three- 
and four-boson systems from our numerical study Q, 
and Fig. [Ijb) shows a schematic representation of the 
region marked by a circle in Fig. [TJa). [Note that we 
use the function F{x) = sgn(a;)ln(l -I- |x|) to visuahze 
the whole energy landscape.] The black solid lines are 
the energies of the four-boson states, the red-solid lines 
represent the dimer-dimer (lower) and dimer-atom-atom 
(upper) collision thresholds, and the green-dashed lines 
are the Efimov trimer energies, or more precisely, the 
atom-trimer collision thresholds. The structure shown in 
Fig. [TJb) repeats every time a increases by the Efimov 
geometric factor e^^^° w 22.7, reflecting the pervasive 
influence of Efimov physics on the four-boson properties. 



Figure [Ijb) illustrates the resonant effects possible in 
dimer-dimer collisions (a > 0). In general, one expects 
pronounced features in ultracold scattering observables 
whenever a bound state crosses, or emerges from, the col- 
lision threshold, or when an additional decay channel be- 
comes energetically available. In Fig. [T]Jb), for instance, 
as a increases and approach a*^, where an Efimov trimer 
is created and causes a divergence in the atom-dimer scat- 
tering length Gad, an infinite number of four-boson states 
Q emerge potentially causing resonant effects in ultra- 
cold atom- molecule gas mixtures Q. For a pure molec- 
ular gas, however, it is primarily collisions of the type 
B2 + B2 that provide the pathway to detailed study of 
four-body physics. As Fig. [TJ^b) shows, there exist pre- 
cisely two four-boson states emerging from the B2 + B2 
threshold at two specific values of a, which we denote 
as a"^^ i (where i ~ \ and 2), causing add to diverge at 
values [l5| . In the vicinity of a^^j, therefore, the dimer- 
dimer interaction can be tuned from attractive, add < 0, 
to repulsive, add > 0, irrespective of the effectively repul- 
sive character (a > 0) of the interatomic interactions. At 
a = a^jj, however, a different process takes place. Here, 
the B2 + B2 and B3 + B channels become degenerate 
[see Fig.[ljb)] and one can expect an enhancement of the 
B2+B2 ^ B3+B rearrangement reaction rate over other 
dimer loss processes. 

A key property resulting from the universality in the 
four-boson problem is that the aforementioned dimer- 
dimer resonances and rearrangement reaction are uni- 
versally related to the three-body Efimov physics. The 
values for a 



dd. 



and a^^ are controlled by a single three- 



body parameter, a*^. We, therefore, express our results 
in terms of the universal ratios 0'*ii4ilo*ad and a'^j^/ a'^j^, 
which determine the positions of the dimer-dimer reso- 
nances, ajjij, and the critical scattering length aj^^ at 
which the B2 + B2 ^ B3 + B rearrangement reaction is 
enhanced. In practice, our calculations extract the values 
for such ratios from the energies in the region near the 
third Efimov state [see Fig. [Ifb)] where the states are 
largely unaffected by nonuniversal short-range physics. 
We obtain 



'(id,! 



2.37, 



^dd,2 



6.6, and -f - sa 6.73. (2) 

"ad "ad '^ad 

We expect that such ratios are universal within an er- 
ror of less than 2%, as determined by comparison of our 
value for a^^/a*^ to the value obtained using the semi- 
analytical results from Ref . Q . 

In Fig. [21 we show our numerical calculations for add 
and for the dimer-dimer relaxation rate, l^^f , calculated 
using: 

k^o Rc[tan5<id] , ,,dd Stt^i ^ 
add = and V.^f = —}_^\Sf. 

f 



-dd 



^(3) 



where fc^ = 2mEcoi/h, with Ecoi = E — 2Eb being the 
collision energy, and exjp{2i6dd) ~ Sdd.dd- The S'-matrix 
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a (a.u.) a (a.u.) 



FIG. 2: (color online), (a) shows poles in the dimer-dimer scattering length add, plotted as a function of the atom-atom 
scattering length a. The poles reflect the emergence of the two universal four-boson states at a = "dd,! ^nd add,2 that are 
associated with every Efimov trimer[see Fig.[TJb)]. These states allow for the control of the dimer-dimer interactions, (b) shows 
the thermally- averaged V^^f and resonant behavior associated with the four-boson states leading to relaxation into deeply bound 
states. For a > add, however, l^^f is dominated hy B2 + B2 B:i + B rearrangement collisions. 



elements are obtained from the solutions of Eq. IT]). In 
our calculations, the only decay channel is the one as- 
sociated with the lowest (relatively deep) Efimov state. 
We, therefore, do not take into account the decay into 
deeply bound two-body states but we still expect 
the general aspects of our results to remain valid. 

Figure [2] shows our numerical results, with the atomic 
mass chosen to be that of Cs, and a*^ = 400 a.u., 
accordingly with the observations of Ref. [§], where 
a*^ = 400 a.u. Figure ^a.) then shows the dimer- 
dimer resonances caused by the four-boson states shown 
in Fig. [ijb). We note, however, that from Fig. [2lja) we 
obtain ct^^/a^^;! ~ 3.57, which deviates from our pre- 
dicted values from Eq. ([2]), o.dd/'^ddi ^ 2.84, in about 
20%. These deviations are produced by finite-range ef- 
fects that introduce nonuniversal corrections to the 
energies of the two-, three- and four-body systems in 
Fig. 1(b). Figure Hfb) shows our results for the thermally 
averaged V^^f ■ clearly exhibits an enhancement of the 



inelastic loss at 



corresponding to the dimer-dimer 



resonance discussed above. Near aj5^2(~ '^dd)^ however, 
the resonance effect is masked by an enhanced rate for 
the rearrangement reaction B2 + B2 — + B3 + B. In 
fact, our calculations show an efficiency of about 98% 
for the B2 + B2 B3 + B reaction. In the pres- 
ence of deeply bound two-body states, we still expect 
a high dimer-trimer conversion efficiency because of the 
stronger overlap between these states than with other 
deeply bound channels. Therefore, the dominance of the 
B2 + B2 ^ B3 + B reaction, and its negligible energy re- 
leased, makes this process promising for the formation of 
Efimov trimers which can (under certain conditions) re- 
main trapped. In a pure molecular sample, a clear exper- 
imental signature of such rearrangement reactions, and 
therefore of trimer formation, would simply be achieved 



by detecting the reappearance of atoms for a > a^^. The 
lifetime of such trimer states can, however, be a critical 
issue, since it depends on its intrinsic decay rate 13, l3| 
and on collisions with other atoms. 

Our results can be compared to the experimental data 
for Cs2 of Ref. [l^l in order to understand the ob- 
served scattering length and temperature dependences. 
In Ref. [13] it was found that V^'^f exhibits a minimum 
for relatively large values of a, and that for the largest 
values of a it was found that ]/^f increases with the tem- 
perature, in contrast to the constant behavior predicted 
by the Wigner threshold law. Based on the temperature 
dependence, which we discuss in details below, we believe 
that the minimum observed in Ref [lo| could be associ- 
ated with our results in the region aj^j ^ < a < a'^^. Even 
though the range of a in which the minimum is observed 
in differs from that observed experimentally, finite 
range corrections [3| can be expected to be more severe 
for the experimental data than for our numerical calcu- 
lations, performed much deeper in the a ^ ro universal 
regime. 

In order to better understand the physics underl yin g 
the intriguing temperature dependence found in Ref. [lO| , 
we present a simple model for relaxation that captures 
the important physics also seen in our numerical calcula- 
tions [Fig. E^b)]. The source of the temperature depen- 
dence in our model, and calculations, is related to the 
presence of a trimer state close to the dimer-dimer colli- 
sion threshold. It is well known that the Wigner thresh- 
old law is applicable only for energies (temperatures) 
smaller than the smallest energy scale in the system [l3l . 
Therefore, any deviation from the Wigner threshold law 
can only be understood in terms of a new energy scale in 
the system that is smaller than the characteristic energy 
scale given by Eh « h/ma? [l^. Based on the knowl- 
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edge of the spectrum for Cs2 (20| , this new energy scale 
can be associated with a trimer state, allowing for the 
enhancement of the B2 + B2 ^ + B reaction as the 
temperature is increased [Fig. Hfb)] . 

In our model, therefore, the important energy scale in 
the system will be A = \E3t — 2Eh\ < Eb, associated with 
the energy of the trimer state. For Ecoi < A, our simple 
model assumes a Wigner threshold law behavior, that 
K-ef ~ hcri/m, where has dimensions of length. For 
Ecoi > A, however, we assume that the rate is unitarity- 
limited, V^if = ^M7r/(25mi)^;7J/^ with < 1. Its 
thermal average is: 



Erf 



A 



A 



2 e-^I^^T" 



hc„ 



(4) 



where ki, is the Boltzman's constant. The inset of Fig. [3] 
shows a comparison between Eq. ([4]) and our numerical 
calculations for a < aj^^ and establishes the validity of our 
analytical model. We then use Eq. ^ to fit the experi- 
mental data shown in Fig.[3]and find excellent agreement. 
However, since our model does not specify the nature of 
the trimer state, we can not determine if the temperature 
dependence in Rcf. is caused by an Efimov state. 

A{a = 800 a.u.) 











A = /r,h - 2Et 
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FIG. 3: (color online). Comparison between experimental 
data (filled circles:' a = 800 a.u.; open circle: a = 500 a.u.) 
and our model [Eq. for the temperature dependence of 
V^f • In our model, the apparent deviation from the Wigner 
threshold law is caused by the presence of a trimer state just 
above the dimer-dimer threshold, as indicated by the vertical 
lines. Inset: Comparison of our model (black-solid line) with 
numerical rates (the red-solid line is our total rate as a func- 
tion of Ecoi, and the green-dashed line shows the thermally 
averaged results as a function of T). 

In summary, we have demonstrated the universal prop- 
erties of ultracold dimcr-dimcr collisions and determined 



their connection with three-body Efimov physics. The 
existence of universal resonances allows for the tunability 
of the dimer-dimer interaction and opens up the possibil- 
ity of study molecular gases with attractive and repulsive 
interactions. Our numerical results show that ultracold 
molecular gases might offer an efficient path for the cre- 
ation of a gas of Efimov trimers via B2 + B2 ^ B3 + B 
rearrangement reactions. In particular, the observation 
of resonant effects in dimer-dimer collisions obeying the 
ratios a'^^^^/a*^ and a'ddlKid^ within the limits imposed 
by nonuniversal corrections we discussed here, is hereby 
shown to provide a simultaneous verification of both 
three- and four-body universal physics. 
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